Semiclassical theory of a quantum pump 
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In a quantum charge pump, the periodic variation of two parameters that affect the phase of the 
electronic wavefunction causes the flow of a direct current. The operating mechanism of a quantum 
pump is based on quantum interference, the phases of interfering amplitudes being modulated by 
the external parameters. In a ballistic quantum dot, there is a minimum time before which quantum 
interference can not occur: the Ehrenfest time te. Here we calculate the current pumped through 
a ballistic quantum dot if te is comparable to the mean dwell time td. Remarkably, we find that 
the pumped current has a component that is not suppressed if te S> td. 

PACS numbers: 73.23.-b,05.45.Mt,73.50.Pz 
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I. INTRODUCTION 

A suitably chosen periodic perturbation of an elec- 
tronic device may result in a direct current, even when 
the system is not biased. In the adiabatic limit, when 
the applied perturbation is slow in comparison to the 
escape rate to external contacts, the electronic state 
is the same after each period while a finite charge is 
transferred through the device during each cycle. Such 
devices are called charge pumps. Charge pumps have 
been proposed as current sources and minimal-noise cur- 
rent standards?!*^ as well as diagnostic tools to moni- 
tor how mesoscopic devices respond to changing external 
parameters J££&2&±2 

In this article, we consider a so-called 'quantum charge 
pump'. In a quantum pump, the external perturbation 
affects the phases of the electron wavefunctions only, not 
the classical dynamics of the electrons ) 4 ' '^' 9 ' ' 12 ' 13 Such 
a scenario can be realized in a ballistic quantum dot, 
where small variations of shape-defining gate voltages or 
of an applied magnetic field are known to significantly af- 
fect the electronic phase while leaving classical trajecto- 
ries unaffected. Such a device, with two shape-distorting 
gate voltages to drive the current, was built by Switkes 
et alA Although the pumped current in the original ex- 
periment of Ref. was obscured by rectification effects^ 
'quantum pumping' was demonstrated in a later experi- 
ment, be it not in the adiabatic regime.— 

Microscopically, the current in a quantum pump is the 
result of quantum interference: two classical trajecto- 
ries initially a microscopic distance apart (a Fermi wave- 
length) split and rejoin, with a phase difference that de- 
pends on the applied perturbation. Indeed, if the mean 
dwell time td in the device exceeds the dephasing time 
T<h, the pumped current is suppressed ) 16 i 17 Aleiner and 
Larkin have pointed out that, in a ballistic quantum dot 
(as in any quantum system with a well-defined classical 
counterpart for which the dynamics is chaotic), there is 
another time that governs quantum interference : 18 ' 19 the 
Ehrenfest time te. The Ehrenfest time is the time dur- 
ing which two classical trajectories initially a Fermi wave- 
length apart are separated to a macroscopic distance (the 
width of the contacts or the size of the quantum dot) by 



the chaotic classical dynamics . 20 ' 21 
is given by 



The Ehrenfest time 



In N + const, 



(1) 



where A is the Lyapunov exponent of the classical dynam- 
ics in the quantum dot and N denotes the total number 
of channels in the contacts. The constant in Eq. JT]) is 
independent of N. While the dephasing time governs 
the suppression of interference for long dwell times, the 
Ehrenfest time is the minimum time needed for the ap- 
pearance of interference effects. 

Being an interference effect, the current pumped 
through a ballistic chaotic quantum dot is a random func- 
tion of any parameter that affects quantum phases, such 
as the shape of the dot or the applied magnetic field 
(taken at a reference point during the cycle). There- 
fore, one typically considers averages and fluctuations of 
the current, taken with respect to an ensemble of dots 
with slightly varying shape or with respect to a range 
of magnetic field strengths. Since the ensemble- averaged 
current (I) = 0, the magnitude of the pumped current 
is measured through the second moment (I 2 )- With a 
few exceptions ) 13 ' 22 random matrix theory has been the 
preferred framework for a statistical theory of quantum- 
dot based charge pumps. 11 ' 23 ' 24 ' 25 ' 26 ' 27 Random matrix 
theory describes the regime te < td, in which quantum 
effects occur on a time scale much shorter than the mean 
dwell time td in the quantum dot. In this article, we are 
interested in a quantum pump in the regime in which te 
and td are comparable. Since, typically, tdA ^> 1, this 
requires that we consider the classical limit N ^> 1. 

A remarkably wide range of TE-dependences has been 
reported in the literature for various quantum inter- 
ference effects in ballistic quantum dots. Weak lo- 
calization, the quantum interference correction to the 
ensemble-averaged conductance in the presence of time- 
reversal symmetry, is suppressed cx exp(— te/td) if te is 
large . 18 ' 28 ' 29 ' 30 On the other hand, conductance fluctua- 
tions have no te dependence ) 31 ' 32 ' 33 whereas the quan- 
tum correction to the spectral form factor of a closed 
quantum dot with broken time-reversal symmetry exists 
for finite te only ; 34 ' 35 being absent in the random matrix 
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limit te — > Oi2& The main result of the present article is 
that current pumped in a quantum pump has yet another 
TE-dependence, 



(I 2 ) = k± + 



-2t e /t d 



W 2 



'RMT, 



(2) 



where (/ 2 )rmt is the prediction of random matrix theory. 

We calculate the Ehrenfest-time dependence of the 
pumped current using a semiclassical theory in which 
the pumped current is expressed as a sum over classical 
trajectories in the quantum dot. Our theory goes be- 
yond a previous semiclassical theory of quantum pumps 
by Martinez-Mares et al.^ who did not consider the 
role of the Ehrenfest time. Whereas Ref. calculated 
the pumped current in the 'diagonal approximation', in 
which only pairs of identical classical trajectories are con- 
sidered, we also include the leading off-diagonal terms in 
the sum over classical trajectories. In this respect, our 
theory builds on previous work by Richter and Sieber 37 
and Heusler et al-r^ who developed a systematic way to 
include off-diagonal terms in trajectory sums of the weak 
localization correction to the dot's conductance and other 
quantum interference corrections to transport. Unlike 
the diagonal approximation, which is known to violate 
current conservation, the present version of the semiclas- 
sical theory is fully current conserving. 

In Sec. |TI] we express the pumped current in terms of 
a sum over classical trajectories. This section closely fol- 
lows Ref. The sum over classical trajectories, which 
gives the final result |(3J), is then performed in Sec. IIIII 
We compare the theoretical predictions to numerical sim- 
ulations in Sec. IIVI and conclude in Sec. [V] 



II. SEMICLASSICAL THEORY OF A 
QUANTUM PUMP 

The geometry under consideration is shown in Fig. [T] 
For definiteness, we consider a ballistic chaotic quantum 
dot coupled to two electron reservoirs via ballistic point 
contacts. The quantum dot has two holes. Slow peri- 
odic variations of the magnetic flux through each of the 
holes change the phases of wavefunctions, not the classi- 
cal trajectories. Both holes are 'macroscopic': their size 
is large in comparison to the Fermi wavelength \p and 
their boundary is smooth. 

We assume that the variation of the fluxes is adiabatic, 
with period 2it [uj much larger than either the dwell time 
or the Ehrenfest time. In the adiabatic regime, the time- 
averaged current II through the left contact can be ex- 
pressed in terms of the N x N scattering matrix S of the 
quantum dotji 1 - 



II = — 



eu> 
2^ 



n 1 



2 n L (4,!,h). 



N Ni 



(3) 
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FIG. 1: Left: Schematic drawing of the quantum pump un- 
der consideration. The quantum pump consists of a ballistic 
quantum dot with two holes, through which a time-dependent 
magnetic flux <f> is applied. The classical dynamics in the dot 
is assumed to be chaotic. Right: pumping cycle in flux-space. 



Here <pj is the flux through hole j, j = 1,2, measured in 
units of the flux quantum hc/e, and the integral in Eq. 
© is taken over the area enclosed in (^q, <^ 2 )-space in one 
pumping cycle, see Fig.[TJ The number of channels in the 
left and right contacts are Aq and N%, respectively, and 
TV = Aq + A^. In Eq. §5§ we assumed zero temperature. 
The current through the right-contact reads 

Ir = ^ / dep^n^uh), 

N N 

n fl (^,0 2 ) = -im]T Yl 



Current conservation implies 

I = Il = -Ir, 

so that 



IT 



TT 



(4) 



(5) 



(6) 



The key point of the semiclassical approach is to ex- 
press the scattering matrix S as a sum over classical tra- 



jectories 7 ) 39 ' 40 i 41 



(J) 



where 7 starts at a contact with transverse momen- 
tum compatible with lead mode v and ends with trans- 
verse momentum compatible with mode /i. (The modes 
fi = 1, . . . , Ni are in the left lead and the modes /i = 
N± + 1, . . . , N are in the right lead.) Further 6> 7 is the 
classical action (which includes Maslov phases) and A 7 is 
the stability amplitude. The classical action is modified 
by the presence of the fluxes ipi and 2 as 



= 5 ( 7 ) + /#iWi( 7 ) + 02W 2 (7)], 



(8) 



where Wj(j) is the winding number of the trajectory 7 
around the hole j, j = 1,2. We'll be interested in the 
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regime <fij <C 1, and Wj ^> 1, for which the discreteness 
of the winding numbers does not play a role. Substituting 
Eqs. © and © into Eqs. (3} and (gj), we arrive al^ 



N Ni 



it 



EE E A 7lJ 4 72 Wi(7i»72) 



4-7T 

JVtd 

f=l jU=l 7l<72 

x sin[(5 72 -S 71 )/?i], 

N N 

4ir 



E E E ^^^1(71)^2(72) 



f=l At=IVi + l 71,72 

x sin[(5 72 -5 71 )/ft], 



(9) 



where both 71 and 72 are compatible with modes ^ and 
fi upon entrance and exit, respectively. 

In the remainder of this article we limit ourselves to 
bilinear response, I oc A. For bilinear response, the 
fluctuations of the kernels H L and 11^ can be neglected 
when performing the integral over the enclosed area A 
in (<pi, 2 )-space, so that it is sufficient to calculate the 
average and variance of the kernels H L and 11^ in order 
to find the average and variance of the pumped current 
I, 



{I) 
(I 2 ) 



etuA 



IF 



ecuA 
~2^T 



4tt 2 



-((n 1 



e 2 u 2 A 2 
4tt 2 



n L u R ) 



(10) 

(11) 



In the semiclassical framework, this means that we can 
neglect the dependence of the classical actions <S 7j . on the 
variations of the parameters <pi and 4>2 in the sines in Eq. 
©■ 

Before we can perform the summation over classical 
trajectories necessary to calculate (I 2 ) we must specify 
the statistics of the winding numbers W\ and W%. For 
a quantum dot with chaotic classical dynamics, one may 
assume that the two winding numbers W± (7) and W% (7) 
are statistically uncorrelated, and that the average wind- 
ing number is zero, 

(Wj(y))=0, j = 1,2. (12) 

For the variance of the winding number we take^ 2 - 

(W j fr) 2 ) = C j t, j = l,2, (13) 

where the coefficients Cj depend on the size and shape 
of the quantum dot and on the location of the hole j, 
j = 1,2, but not on the Ehrenfest time and t is the du- 
ration of the trajectory 7. For a chaotic quantum dot, 
one may also assume that the winding numbers Wj(7i) 
and Wj(72) of different trajectories 71 and 72 are un- 
correlated, j = 1,2, except if there exist strong classical 
correlations between the two trajectories. Strong classi- 
cal correlations exist if the two trajectories are within a 
phase space distance c, where c is a cut-off below which 
the chaotic classical dynamics in the quantum dot can 
be linearized. If two trajectories are correlated for only a 




FIG. 2: (Color online) Schematic drawing of two partially 
correlated trajectories. The true trajectories are piecewise 
straight with specular reflection off the dot's boundaries. 



part of their duration, as shown schematically in Fig. [21 
only the duration of the correlated segment contributes 
to the average of the product of the winding numbers, 



(W j (j 1 )W j ( l2 ))=C j t b , 3 = 1,2. 



(14) 



Although our theory is formulated for a quantum 
pump where the current is driven by time-dependent 
magnetic fluxes, the mathematical framework to describe 
other geometries, e.g. a pump in which current is driven 
by shape changes, is identical. 22 For generic perturba- 
tions of a chaotic quantum dot the parameters still en- 
ter the classical action as a Gaussian random process, 
and the actions of different trajectories will have statis- 
tically uncorrelated dependences on the parameters. In 
that sense, our theory applies to arbitrary perturbations, 
not only time-dependent magnetic fluxes. (The only ex- 
ception is a pump in which a gate voltage is changed 
uniformly in the quantum dot. A similar separation be- 
tween 'generic' perturbations and a uniform gate voltage 
shift exists in the random matrix description of quantum 
pumps. 



1 i :1 1,26-) 



III. SUMMATION OVER CLASSICAL 
TRAJECTORIES 

The semiclassical framework and the statistics of the 
winding numbers outlined in the previous section is the 
same as that used in a previous semiclassical theory 
of quantum pumps by Martinez-Mares, Lewenkopf, and 
Muccioloi^ 2 - In Ref.[23, the fourfold sum over classical tra- 
jectories required to calculate the current variance is eval- 
uated in the 'diagonal approximation', in which the four 
trajectories contributing (I 2 ) are taken pairwise equal. 
Although Ref. reports agreement between semiclas- 
sics and random matrix theory, the use of the diagonal 
approximation to calculate the current variance (I 2 ) is 
problematic because it violates the unitarity relation ([6]) : 
In the diagonal approximation, one has (II L II fl ) = 0, 
while ((II L ) 2 ) ^ 0. The diagonal approximation also fails 
to account for the fact that the variance of the pumped 
current does not depend on the presence or absence of 
time-reversal symmetr y 11 ^ 3 

Problems with the diagonal approximation are not lim- 
ited to the semiclassical theory of a quantum pump. In 
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FIG. 3: (Color online) Schematic drawing of trajectory con- 
figurations of types A (top) and B (bottom) that generate the 
pumped current. Type A consists of all quadruples of classical 
trajectories with two separated small-angle encounters. Type 
B consist of all quadruples of classical trajectories for which 
two (or more) of the trajectories revolve around a periodic 
orbit. If the Ehrenfest time is much smaller than the mean 
dwell time, the pumped current is carried by trajectories of 
type A. If te S> to, the current is carried by trajectories of 
type B. 



fact, the diagonal approximation is well known to fail 
at providing a current-conserving description of other 
quantum interference effects, such as the weak localiza- 
tion correction to the conductance and universal con- 
ductance fluctuations 4 s . For weak localization, this prob- 
lem was solved by Richter and Sieber, who were able to 
include the relevant off-diagonal configurations of clas- 
sical trajectories into the trajectory sum^i Below, we 
use a generalization of Richter and Sieber's technical 
innovatio n 33 ' 38 ^ 4 ' 45 to perform the trajectory sums re- 
quired for a current-conserving semiclassical theory of a 
quantum pump. 

Since the two winding numbers W\ and W2 are statis- 
tically uncorrelated and have zero average, the trajectory 
sum in Eq. ^ immediately gives (H L ) — (H R ) = 0, and 
hence 



(/) = 0. (15) 



In order to calculate the average of the squared current, it 
is technically most convenient to use the second equality 
in Eq. (JTT|) . Using the semiclassical expression one 




FIG. 4: (Color online) Trajectories of type A without and 
with correlated entry into the quantum dot (top and bottom, 
respectively) . 

has 

2 N JVi N 

D j_L,p—l a— JVi + 1 71,72)73)74 
x \A lx Aj 2 Asy 3 Ary 4 

x Wi (71)^(72) sin[(S 72 - 5 7l )/fi] 

x (73)^2(74) sin[(5 74 - SyJ/h]) .(16) 

The non-vanishing contributions to the current variance 
are from quadruples of trajectories for which the two ac- 
tion differences in the sine functions are almost identical, 
resulting in a non-oscillatory contribution to the trajec- 
tory sum. However, this alone does not ensure a finite 
contribution to (I 2 ): In addition, the ensemble averages 
of the winding numbers should not vanish. This means 
that the trajectory pairs 71 and 73, and 72 and 74 should 
be correlated for (at least) part of their duration. How- 
ever, since these trajectories exit the dot in different con- 
tacts, they can not be identical. 

The configurations of classical trajectories that meet 
these requirements are the same as the trajectories 
that contribute to conductance fluctuations^ 3 - They are 
shown in Fig. [3] The first type (labeled "A" ) consists 
of two pairs of trajectories that undergo two separate 
small-angle encounters. Between the encounters and the 
contacts 71 and 72, and 73 and 74 are pairwise equal 
(up to quantum uncertainties) . Between the two encoun- 
ters, 71 is paired with 73 or 74, while 72 is paired with 
the remaining trajectory, 74 or 73. The second type of 
trajectories (labeled "B") contains a closed loop, which 
is traversed one (extra) time by two of the four trajec- 
tories. Below we discuss both contributions separately. 
Since the calculations closely follow the calculations of 
the conductance fluctuations, we refer to Ref. l33| for de- 
tails of the formalism, and here restrict ourselves to those 
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parts of the calculation where we differ from Ref. l33l . 

A. contributions of type A 

Within the configurations of type A, one distinguishes 
trajectories for which the first small-angle encounter fully 
resides inside the quantum dot, such that the entrance of 
the pairs trajectories 71 and 72, and 73 and 74 is uncor- 
rected, and a configuration of trajectories in which the 
first small-angle encounter touches the entrance contact, 
so that all four trajectories enter the dot together. The 



two situations are shown in Fig. [4j Since 71 and 72 exit 
through the left contact, whereas 73 and 74 exit through 
the right contact, the second encounter must fully reside 
inside the quantum dot in all cases. We refer to the two 
cases shown in Fig. [4] as "Al" and "A2". 

The durations of the encounters, defined as the time 
during which the phase space distance between the tra- 
jectories is less than a classical cut-off scale c, are denoted 
i nc,i and i e nc,2, as depicted in Fig. [4] The durations of 
the uncorrelated stretches between the encounters are de- 
noted by n and T2 . We then expand 



4sin[(S 72 - S 71 )/h] sin[(5 74 - S l3 )/h] = e^-S^+S^-S^/h _ e i(s 72 -s 71 +S 74 -S 73 )/ft + c c ^ 

and find quadruples of trajectories 71, 72, 73, and 74 for which the exponents in Eq. (|17[) are small. Between the 
encounter regions and the contacts, 71 and 72, and 73 and 74 are paired, resulting in only a small contribution to 
the action differences in Eq. (fTT|) . In order to ensure that the contribution from the stretch between the encounters 
is small as well, one pairs 71 with 73 and 72 with 74 between the encounters for the first term in Eq. (|17p . and 
71 with 74 and 72 with 73 for the second term. In the former case {W\ (7i)Wi (73)) = Ci(i onCi i + i nc,2 + Tl) and 
(W 2 (72) W 2 (74)) = C 2 (Wc,i + t C nc,2 +t 2 ). In the latter case, one has (Wi (71)1^1(73)) = Ci(i onc ,i + t cnc , 2 ) and 
(W2 (72) 1^2(74)) = C 2 (t enCt i + icnc.2)- At each of the encounters we take a Poincare surface of section. The action 
differences AS = <5> 72 — 5 7l + 6> 73 — 5 74 or AS — 5 72 — 6> 7l + 5 73 — <S> 74 that appear in the exponents in Eq. (jTTJ) 
are then expressed in terms of the phase space distances Sj and Uj between the two solid trajectories in Fig. [3] along 
the stable and unstable directions in phase space at each of the encounters, AS — siUi + s 2 u 2 i 46 ^ 47 The encounter 
durations < C nc,i and i cnCj 2 are expressed in terms of these coordinates as 



*encj = J^T^Tl, J = 1,2. (18) 



Calculating the trajectory sum as in 

Refs.mmsnE 

we then find 



enc,l^cnc,2 

xC*iC 2 [(tenca + ic„c :2 )(n + r 2 ) + t x t 2 \ . (19) 

The factors £ cnc ,i and i on c,2 in the denominator in Eq. (|19p cancel a spurious contribution from the freedom to take 
the Poincare surface of section at an arbitrary point in the encounter region. 

We then make the variable change Uj — c/<jj, Sj — exjdj. In terms of the new variables, one has t onc j = 
A -1 ln(l/|xj|). Performing the integrals over t\, t 2 , a\, and a 2 , one then finds 



„' 1 «0 



2 r rl 

[27b (i enc, 1 + ^onc,2) + 7"d] , (20) 

J'=l ' 

where we abbreviated r — c 2 /%. 

We now turn to compute contributions of type A2. Because the first encounter touches the leads, the encounter 
duration i cnc ,i is no longer determined by the phase space coordinates s\ and u\, but becomes an integration variable 
itself. The integration range for t en c,i is A -1 ln(c/|ui|) < i cnc ,i < A -1 ln(c 2 /|uiSi |), 45 so that 




e -l(sini+S 2 M2)/'i-(tonc,l+ t onc,2)/TD 



X ■ 



{2irh) 2 t enCil t enCt2 



-C X C 2 [(t enc ,i + t onCi2 )(ri + r 2 ) + t x t 2 \ . (21) 
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The second encounter time is still given by Eq. (fTS|) above. As before, we perform the variable change Uj = c/<?j, 
Sj = cXjGj, j = 1,2. Upon integration over t\, t 2 , a%, 02, and i cn c,i we then find 



N 2 



rroA / dx2X 2 TD cos(rx2) 



x rroA / c?xi cos(ra;i) 



1/Atd 



(2r D i c 



2td£ 



3rg)-(2r D t c 



(22) 



Here we reverted to the notation i cnc ,i = A -1 ln(l/|xi|) in order to make contact with the result for (H l H r )ai 
obtained previously. The second term between the square brackets [. . .] in the second line of Eq. (j22|) leads to a 
rapidly oscillating function of r and will be neglected. Combining the remaining terms with Eq. (|20[) and taking the 
limit h — > (r — > oo) while keeping the ratio te/td fixed, we find 



1 -|2 

rA / dxx 1 ^ 7 " cos{rx) 



ifi 2 ^1^2 ^ ~ 2 -ja 
~ 167r -p-^iCarge t d . 



(23) 



where the Ehrcnfcst time is defined as 

^ = Ilnr=iln^. (24) 

Note that the limit 7i — * at fixed te/td, taken above, 
can be realized by narrowing the lead openings such 
that the dwell time increases. While this procedure may 
change the trajectories contributing to pumping, their 
statistical properties should be (almost) unaffected. In 
particular, the Lyapunov exponent will converge to its 
value in the closed cavity. 

B. contributions of type B 

The trajectory configurations of type B provide the 
dominant contributions to the pumped current at large 
Ehrenfest times. Each pair (71,72) and (73,74) consists 
of a 'short' trajectory and a 'long' trajectory, where the 
long trajectory differs from the short trajectory by wind- 
ing once around a periodic orbit, see Fig. [3J (Strictly 
speaking, one should say that the long trajectory winds 
one extra time around the periodic orbit, because all tra- 
jectories involved will wind multiple times around the 
periodic orbit if te is much longer than the period t p 
of the periodic orbit.) In the trajectory sum (|9|), e ach 
trajectory can be the short one. Following Ref. l33l. we 
parametrize these configurations by Poincarc surfaces 
of sections which measure the phase space coordinates 
(si, Ui) and (S2, 112) of two short trajectories with respect 
to the stable and unstable manifolds of the periodic orbit. 

It will be useful to define several relevant times. The 
two short trajectories in the configuration are depicted 
in Fig. [5l The time for which the short trajectories are 
correlated with the periodic trajectory (phase space dis- 
tances smaller than the cutoff c) is denoted £ C nc,i [for the 
short trajectory of the pair (71,72)] and t nc.2 [for the 
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FIG. 5: (Color online) Definitions of the relevant durations 
for trajectories of type B. The figure shows the short trajec- 
tories from each pair, as well as the periodic orbit. The long 
trajectories from each pair, which are not shown, differ from 
the short trajectories by one revolution around the periodic 
orbit. The figure shows the durations i e ac,i and t C nc,2 of the 
encounters of the short trajectories with the periodic orbit, as 
well as the time t s that the trajectories are correlated before 
arriving at the periodic orbit. A similar time r u measures 
the time the trajectories are correlated after departure from 
the periodic orbit. In the example shown in the figure, the 
trajectories depart at different points, so that t u — 0. 



short trajectory of the pair (73,74)]- By our construc- 
tion, the longer trajectory in the pair (not drawn in Fig. 
[5]) will be correlated with the periodic trajectory for a 
time t cnc j + t p , j = 1,2, where t p is the period of the 
periodic orbit. Since the trajectories may leave the neigh- 
borhood of the periodic trajectory together, the correla- 
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FIG. 6: (Color online) The overlap time is denned as the 
overlap of the two encounter times t e nc,i and t C nc,2- 

tions may extend away from the periodic trajectory. If 
applicable, we denote the time for which the trajectories 
are correlated with each other before they arrive to the 
periodic orbit by t s . Similarly, t u denotes the eventual 
correlation time after the trajectories leave the periodic 
orbit. If t s 0, the encounter region extends away from 
the periodic orbit and may touch the entrance contact. 
Hence, we distinguish contributions of type Bl and type 
B2, where type Bl refers to those trajectories for which 
the encounters fully reside in the quantum dot, and type 
B2 refers to configurations of trajectories with correlated 
entry (encounter that touches the entrance lead open- 



ing). Since the trajectory pairs (71,72) and (73,74) exit 
through different contacts, one does not have to con- 
sider the possibility that the encounter region touches 
the exit lead opening. The Poincare surfaces of section 
are taken at points where the short trajectories are cor- 
related with the periodic orbit. The distance between 
the two Poincare surfaces of section is measured by the 
travel time t l2 along the periodic orbit, where we require 

\t 12 \<T p /2. 

Even the shortest trajectory in a pair can wind around 
the periodic trajectory several times. Thus, it is possible 
that t cnc j > T p . In order to separate out full revolutions, 
we write 

where rij is a non-negative integer and < t Gncj < t p . 
The use of t enc j complicates some of the calculations 
since this time is not a continuous function of the surface 
of section coordinates. On the other hand, this time ap- 
pears naturally due to the random Gaussian character of 
the winding numbers W\^ 2 . 

Without loss of generality (but with inclusion of a com- 
binatorial factor two), we take 71 to be one of the short 
trajectories. Expanding the product of sines as in Eq. 
(fP7"|) . we then see that smallness of the the action dif- 
ference requires 73 to be the other short trajectory for 
the first term in Eq. (fT7|) . whereas 74 is the other short 
trajectory in the second term in Eq. (|17|) . In the former 
case, one has 



{Wi (71) Wi (73)) = Ci[nin 2 T p + nit enCt 2 + n 2 W,i + St + t s + t u ] 

(1^2(72)1^2(74)) = C 2 [(m + l)(n 2 + 1)t p + (m + 1)W,2 + (n 2 + l)i C nc,i + Si+ r s + t u }. 

The time St is the overlap of the 'remainders' t enc ,i and i C nc,2, see Fig. El In the latter case, one has 

(Wi(7i)Wi(7 a )) = Ci[m(n 2 + 1)t p + ni i cnc , 2 + (n 2 + l)t „c,i + St + t s + t u ] 
(^2(72)^2(74)) = C 2 [(ni + l)n 2 T p + (m + l)t e nc,2 + n 2 t cncA + St + t s + t u \. 

The resulting contribution of type Bl then reads 

,^Vi^2Td^^ f dr v _ T , f Tp/2 , f c , , , , e^ v ^- U2S ^/ n - { - T ' +T ^/ T ^ 



(Tl L n R ) m = 167r 2 iVl ^f D C 1 C2 / ^e'^ f " dt 12 f d Sl d Ul d S2 du 2 

^ J I'D J-t p /2 J-c 



(2nh) 2 t enCt it enCt2 



X [r p (St + T S + T u ) - tone, 1*0110,2] ■ (26) 

The contributions of type B2 are calculated in a way very similar to the one used for type A2, 

>NiN2T D ^^ f dr p _ T /_ f T » /2 , f c , , , , e^ UlSl - U2S2 ^ h - T "/ TD 



(U L n R ) B2 = l^^l^ Cl C 2 f d ^Le-^™ f dt 12 f d Sl d Ul ds 2 du 2 - 

X { (l - e- r '/ TD ) [ Tp (St + T D + r„) - tonc,i? c „c,2] - e"^r p r s | . (27) 
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Adding both contributions, we then obtain 



»/2 



dt 12 



Tp/2 



dsiduids2du2- 



,i(uisi — U2S 2 )/ft — r u /r D 



{2irh)H 



enc,l^cnc,2 



| U - e Ts / rD ^ T p T D + T p (<5£ + T u ) ~ ?onc,l^onc,2| 

I 



(28) 



In App. [X] it is shown that the contribution from the 
combination r p St — i e nc,iienc,2 between the curly brackets 
{. . .} is of order CiC27d/A, which can be neglected with 
respect to (n^II^)^. The contribution from the term 
t p t u is a rapidly oscillating function of r = c 2 /% and is 
omitted. The only nonzero contribution arises from the 
term t p td, which is the only term that is multiplied by 
both exp(— t s /td) and exp(— t m /td). This factor t p td 
can be taken out of the integral over the phase space 
variables. The resulting integral was computed in Ref. 
l33l . with the result 



(n L n fi ) B = -8tt 

x(l 



,N X N 2 



N 2 



C\C 2 



1 p * 



-2t e /t d n 



(29) 



In Eq. (|29p we omitted terms proportional to e~ XTp . 
Typical periodic trajectories contributing to Eq. (|29|) 
have duration ~ td. In the limiting procedure discussed 
below Eq. (pH)) these contributions scale as exp(— Atd) — » 
0, justifying this approximation. The integral over the 
period of the periodic trajectory is easily calculated, re- 
sulting in 

(I^ri^ = -S^^C^rg (l - e -*W™) . (30) 

Combining both contributions together, we find that 
the total variance of the pumped current is 



N 2 



-A 2 e z uj z CiC 2 T^ 1 



O ' Mi 

e t d 



(31) 



Equation (|3 1 [) agrees with random matrix theory in the 
limit te — > i 11 i 23 i 24 At finite te, the variance of the 
pumped current is reduced below the random matrix 
value, but the reduction is not more than a factor two. 

When te <C td the pumped current is dominated by 
trajectories of type A. For te 3> td the pumped current 
is carried by trajectories of type B. Since such trajectory 
configurations involve a closed loop they are associated 
with fluctuations of the density of states. Although this 
scenario is very similar to that of the Ehrenfest-time de- 
pendence of the conductance fluctuations, the pumped 
current has a TE-dependent part, whereas the conduc- 
tance fluctuations in a chaotic quantum dot are fully te 
independent. The difference occurs, because the mean 
total duration of the trajectories involved in the internal 
loop for trajectories of type A is twice the mean duration 
of the internal loop for trajectories of type B. Conduc- 
tance fluctuations are insensitive to the loop duration (as 



long as the trajectories stay inside the quantum dot), but 
the pumped current is not. 

Equation (131|) is valid for systems with and without 
time reversal symmetry. This is well known in terms of 
random matrix theor y 11 ! 23 We have verified this explic- 
itly in the semiclassical approach, see App. [BJ We also 
verified that calculation of the current variance using the 
correlator ((n L ) 2 ), while technically more involved, gives 
the same result for the pumped current, as required by 
unitarity. 



IV. NUMERICAL SIMULATION 

Despite the remarkable similarity of the semiclassical 
calculation of the pumped current and the semiclassical 
calculation of the conductance fluctuations, the former 
shows a dependence on the Ehrcnfcst time, whereas the 
latter does not. In this section we report numerical sim- 
ulations of the pumped current and compare these with 
our theoretical predictions. 

Because the computational cost of numerical sim- 
ulation of two-dimensional quantum dots with large 
Ehrenfest times is prohibitive, Jacqoud, Schomerus, and 
Beenakker suggested to simulate one dimensional chaotic 
maps instead.— The phenomenology of chaotic maps is 
identical to that of chaotic cavities ) 49 ' 50 but the compu- 
tational cost of simulating a map is significantly lower 
than that of simulating a cavity. Chaotic maps have 
been successfully used to study the Ehrenfest time de- 
pendence of a range of properties of chaotic quantum 

j t„ 29.32.48.51.52.53.54.55 

The quantum map propagates a finite state vector of 
size M = \/2tt% in time, 



^(t + l)=^(t), 



(32) 



where T is the Floquet operator of the map. Assigning 
two consecutive sets of N\ and N2 elements of the vec- 
tor ip to contacts, the map can be used to construct a 
(Ni + 7V 2 )-dimensional scattering matrix S according to 
the rulei&SS 



S(e) =V[l- e ie TQ] 1 e^TV 1 



(33) 



where V is a (N\ + N2) x M matrix projecting on the 
lead sites, Q = 1 — V T V, and e is the quasi energy. 

The map used in our simulations is the 'three-kick 
quantum rotator', 52 which is specified by the Floquet op- 
erator, 



J rt 



(XUY*UYUX) r 



(34) 
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where 



Y 

Xrr. 



i(jM/6tt) cos(2irm/M) 



r -i(M/12ir)V(2irm/M) 
Or, ' 



U mn = M - l ' 2 e-"' A e^ m -^ 2 ' M . (35) 

In this model M is even, but not a multiple of 3. The 
kick potential is given by 

V{9)=K cos(7rg/2) cos 9 + sin(7rq/2) sin 29, (36) 

where q breaks the parity symmetry of the models The 
parameter 7 plays the role of a magnetic field resulting in 
the breaking of time reversal symmetry. This model was 
used in Refs. andHH to compute the weak localization 
correction of the conductance. 

There are two parameters in this model: The chaotic- 
ity parameter K and the 'magnetic field' 7. We will use 
K and 7 as the time-dependent parameters that drive the 
quantum pump. Although both K and 7 appear in the 
classical map and thus affect classical trajectories, the 
variations used in our simulations are sufficiently small 
that the changes of the classical dynamics can be ne- 
glected within the mean dwell time td = M/(N\ + N2). 

In the simulations we took equal channel numbers in 
both contacts, Ni = N2 = N/2. The semiclassical limit 
is taken by increasing the dimension M = \/2itTi of the 
system, while keeping the dwell time td = M/N fixed. 
This means that N is scaled proportional to M. This way 
we are certain that there are no variations in the classical 
dynamics while decreasing %. With these definitions the 
Ehrenfcst time is estimated to be 



Te 



A 



In AT, 



(37) 



plus an iV-independent constant. The Lyapunov expo- 
nent A is calculated independently from a simulation of 
the classical counterpart of the map. 

The pumping strength depends on the two parameters 
Ck and C 7 that describe the rate at which variations of 
K and 7 affect transport. (In the semiclassical theory, 
these parameters were called C\ and Ca.) In order to cal- 
culate Ck and C 7 , we first performed separate numerical 
simulations of the quantities 



VK 



II- 



Ni N 

E E 

o=l /3=iVi+l 
JVi N 

E E 

o=l /3=JVi+l 



dS a(: 



dK 



dS, 



a 11 



d-, 



(38) 



and compared these to the result of a semiclassical cal- 
culation, 



(39) 





FIG. 7: (Color online) The normalized current variance R, 
as a function of N. The symbols denote simulation results 
for the three-kick model with K — 10, q = 0.2 and 7 = 5 
(circles) or 7 = f (squares). Results for td = 5 (left panel) 
and td = 10 (right panel) are shown. The solid and dashed 
line are fits to theory, as described in text. 



We then calculated the kernels H L and numerically, 
and considered the ratio 



R 



Nj_n 2 /n L n fl 



8tt 2 \ y K y~i 
The prediction of the semiclassical theory is 



R=l 



3 -2t e /t d 



(40) 



(41) 



The normalized pumping current variance R is ensem- 
ble averaged by variation of the quasi energy e and lead 
positions. Since the variation of the lead positions may 
result in large classical variations, we use the same set 
of lead positions for different values of M. The factors 
n L n fi , dk and y 7 are first averaged over the quasienergy. 
This average can be done analytically.™ Then R is aver- 
aged over 2 000 different lead positions. 

We have checked that when the dwell time is increased 
the pumping current variance tends to the random ma- 
trix limit R = 2. (Data not shown. The trend can be 
seen by comparing the left panel of Fig. [7] [where td = 5] 
to the right panel [where td = 10].) However, in order 
to access the regime in which te and td are comparable, 
we have to take the rather small value of td = 5. In this 
case, full ergodicity is not achieved during the time td, 
which results in a deviation from the random matrix limit 
even for small Ehrenfest times. A similar deviation from 
random matrix theory was observed in numerical sim- 
ulations of weak localization and universal conductance 
fluctuations^^^^^^^ 

While the necessity to consider short dwell times pre- 
vents a quantitative comparison of simulations and the- 
ory, our simulations, which are shown in Fig. show an 
unambiguous dependence on rg. To quantify the depen- 
dence we have fitted the simulation results to a function 
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of the form R = A + B exp[-(2/Ar D ) In AT]. The values 
of the Lyapunov exponent are obtained from a simula- 
tion of the classical map££ corresponding to ([34]) . Wc 
find A ~ 1.43 for 7 = 5 and A ~ 2.04 for 7 = 10. The 
fit constants are found to be A ~ 0.74 and B ~ 1.18 for 
7 = 5 and r D = 5; A ~ 0.64 and B ~ 1.0 for 7 = 10 and 
td = 5 ;A ~ 0.56 and B ~ 1.37 for 7 = 5 and td = 10; 
and A ~ 0.61 and B ~ 1.15 for 7 = 10 and r D = 10. 
Deviations from the theoretical values A = B = 1 are at- 
tributed to the nonuniversal corrections (caused by the 
small dwell time To). The unknown additive constant in 
the expression (f3"T)) for the Ehrenfest time is absorbed in 
the parameter B and may cause further deviations from 
the predicted value 5 = 1. 

As mentioned above, the small dwell times used in the 
simulation cause non-universal and systematic deviations 
from the theoretical prediction. Increasing the dwell time 
is numerically expensive. It also strongly suppresses the 
dependence of R on the number of channels N. With 
these limitations, we conclude that, while the simulation 
does not unambiguously confirm the semiclassical theory, 
it certainly is not inconsistent with it: The data shows 
a clear dependence on the Ehrenfest time, whereas the 
nonzero fitted values for A clearly point to a finite value 
of R in the classical limit N — > 00. 



V. CONCLUSION 

In this article we presented a semiclassical theory of 
the current pumped through a ballistic and chaotic quan- 
tum dot with two time-dependent external perturbations. 
In a quantum pump, the external perturbations couple 
to the phase of the electronic wavefunction, not to the 
classical dynamics, so that the pumped current is gener- 
ated through quantum interference only. Our semiclas- 
sical theory identifies the two classes of interfering tra- 
jectories that contribute to the root-mean-square current 
rxasl — {I 2 ) 1 / 2 . These are shown in Fig. [31 In each case, 
the semiclassical contribution for (I 2 ) originates from two 
pairs of classical trajectories which meet at small-angle 
encounters. The interference arises from time-dependent 
phase differences accumulated between the encounters. 

Our semiclassical calculation shows that the action of 
the external perturbations during the small-angle en- 
counters does not generate a direct current. This fol- 
lows, e.g., from the final expressions (f2"3")l and (|29p for 
the two contributions to the pumped current, which con- 
tain the mean times of the interference loops but not the 
encounter times. Although such a fact appears natural 
for a quantum interference effect — nontrivial interfer- 
ence can occur only if trajectories are separated in phase 
space — , within the semiclassical formalism it follows 



only after cancellation between contributions of trajecto- 
ries with and without correlated entry. In this respect, 
our calculation again emphasizes the necessity to treat 
trajectory configurations with and without correlated en- 
try on the same footingi2£i2i^i£& 

The main goal of our calculation was to provide a quan- 
titative description of the pumped current if the Ehren- 
fest time te is comparable to the mean dwell time td 
in the quantum dot. The Ehrenfest time is the time re- 
quired for the chaotic classical dynamics to separate two 
phase space points initially a Fermi wavelength apart to 
classical distance. As such, te presents a threshold time 
before which no quantum interference can occur. While 
we do find that the pumped current is sensitive to the 
ratio te/td, the pumped current remains nonzero, even 
in the limit te S> td. In the limit te 3> ttj, the trajec- 
tories that contribute to the pumped current wind many 
times around a periodic orbit. The escape probability for 
such a configuration depends on the period of the peri- 
odic orbit, not on the total duration of the trajectories 
involved. This explains why such trajectories continue 
to exist even if their duration far exceeds the mean dwell 
time td- 
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APPENDIX A: REMAINDER TIMES 
INTEGRATION 

In this appendix we show that the integral 

f - 1 f Tp/2 f 
J = I dr p e Tp ' TD / dt\2 \ ds\du\ds2du2 

J J-t p /2 J-c 

e i(u 1 S 1 -U 2 S 2 )/h-T s /T D 

X -5 \jp5t - ienc,l*enc,2] >(A1) 

'£ ^enc,l^enc,2 

is of order td /A. After the replacement of r s by t u , which 
is allowed by symmetry, this proves the claim that the 
terms proportional to T p 5t— i e nc,i£enc,2 in the curly brack- 
ets {. . .} in Eq. (|28p does not contribute to the correlator 

(n L u R ). 

First, we point out that the integral over ti2 vanishes 
if one omits the factor exp[— t s /td] from the expression 
for J. Hence, we may replace the factor exp[— t s /td] in 
Eq. (TMJ by exp[-r s /r D ] - 1, 
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/l-Tp/2 re e i(uis 1 -u 2 s 2 )/h . . 
dT p e- T " /TD dt 12 ds l du 1 ds 2 du2— i [e- T ^ TD -l ) (r p <Jf- i cnc ,if onc , 2 ) . (A2) 

The contributions for this integral are from trajectories 71, 72, 73, and 74 with correlated arrival at the periodic 
orbit, since the integrand vanishes if r s = 0. On the other hand, the departure from the periodic orbit need not be 
correlated. 

Following Ref. [H we perform a change of variables using Sj — c/aj, Uj — cXjOj, j = 1,2, and t' 12 = t\ 2 — 
A -1 In (c/|si|) + A -1 In (c/ 1 S2 1 ) • This change of variables amounts to choosing the Poincare surfaces of section at the 
point where each trajectory first comes within a phase space distance c from the periodic orbit. The time i' 12 is the 
distance between the new Poincare surfaces of section, measured as the travel time along the periodic orbit. The 
integrations over aj, j = 1, 2, can be done, and cancel the factors t cnc j in the denominator. The resulting expression 
for J is 




where the sum over the signs originates from the possi- 
ble relative signs between S\ and s 2 and r — c 2 /h. We 
perform another change of variables^ defining u\ = x%, 
u 2 = ±x 2 and w = ±e A * 12 . Both the integration over t' l2 
and the summation over the signs are combined into the 
integral over w. A straightforward calculation leads to 



J 



-Is. f dw 
ar„e T o / - — - 
w\ 



du\du 2 



-t,/t d 



1 



x cos (rui - ru 2 ) (r p 8i - t cnC7 it cnc , 2 ) 



(A4) 



In Eq. (|A4p . the time t s depends only on the loca- 
tion where the trajectories in each pair become correlated 
with the periodic orbit. Following Ref. [H, we set 



In 



max(u>, 1/w) — 1 



6-1 



(A5) 



where b is a numerical constant of order unity. Because of 
the factor exp(— t s /td) — 1, only \w\ close to unity, corre- 
sponding to A|i' 12 | < 1, contributes to the integral. This 
implies that the overlap time St can be approximated by 



St = min(i 



cnc.l ) ^0110,2 ) : 



where 



oncj = -ln(l/|Wj|) 



mod 



(A6) 



(A7) 



The terms neglected in this approximation can, at most, 
lead to corrections of order l/Ar D . 

The integrals over w and U\, u 2 factorize. Integrat- 
ing over w gives a factor of order 1/Atd- The remaining 
integral over u\ and u 2 is complicated due to the dis- 
continuity of the encounter times when they reach the 
period t p , see Eq. (|A7|) . We deal with this problem by 
separating the integral into domains where the integrand 
is continuous. The integral of the term proportional to 
t e nc i^enc,2 factorizes to a product of two integrals. They 



are given by 



r j duj cos(rUj)t cnc j 



1 f 1 diij 



sm(r Uj ) - r p /(r), (A8) 



where 



f(r) 



OO 

E 

m=0 



;in (r 



(A9) 



In the classical limit r — > 00 and the integral in the first 
term in Eq. (|A8|) converges to ir/2. This leads to 



■h. 



n 
2A 



r P f(r) 



(A10) 



For typical periodic orbits Ar p 3> 1 and thus we expect 
the second term in (|A10[) to dominate. Since both r and 
e Ar p are j ar g e we expect that the terms in the series /(r) 
are either rapidly oscillating or small. Averaging over 
variations of the cut-off r or a local average over A or r p 
will smear out the oscillating terms. In this case the series 
converges to some small value. We can conclude that / 
is at most of order unity, so that J±j is at most of order 
t p . This suggests that the product Ji,i Ji,2 contributes 
a term of order r p to J. However, we will find that 
the terms proportional to r p cancel and the integral is 
actually smaller. 

It remains to calculate the integral of the term propor- 
tional to TpSt in Eq. (|A4[> . 



J 2 =r\r 2 j du\du 2 cos{r\U\) cos{r 2 u 2 )T p 5t. (All) 

Again, we separate the integration domain into regions 
where the integrand is continuous. Rescaling the integra- 
tion variables as Uj = yje~ mXTp , one finds 
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r\r 2 
A 



V r p e-(" l +')^ / dyi / dy 2 cos (ne^'ji) cos (r 2 e- ZA ^y 2 ) ln yi 



i,m=0 



( dy x ( dy 2 cos(r ie mXTp yi) cos (r 2 e ZATp y 2 )ln?/ 2 

J e - XT P Jyi 



(A12) 



Here we used two a priori different cut-offs n and r 2 for the two integrations, instead of the single cut-off r = c 2 /%. 
The use of two different cut-offs is needed to exclude products of oscillating functions of n and r 2 , which would give 
a non-oscillating contribution to the integral if the same cut-off would have been used. The integral over y 2 in the 
first term of (|Xl2| can be easily done. Similarly, the integral over y\ can be done in the second term, with the result 



dy cos (r ie - mXT "y) sin (r 2 e- tXT "y) - sin [r 2 e- (l+1)XT ^ 



A 



Z" 1 r / 

T p e- lXT " / dy cos {r 2 e- lXT "y) sin (r ie - mXT "y) - sin r ie - (m+1)AT ' 



,i=0 



In y 
lny. 



(A13) 



Partial integration results in 



J2 = E T W i sin (rie-^y) sin (rae"'^!/) + r p 2 /M/(r 2 ) - [/(r x ) + /(r 2 )] , (A14) 



l,m=0 



where the function f(r) was defined in Eq. (|A9|) above. 

In order to estimate the series of integrals in Eq. (|A14p it is useful to divide the double sum into parts using 



oo oo 



oo oo 



l,m=0 l=m=0 m=0l=m+l 1=0 m=l+l 

and consider each part separately. The diagonal sum gives 

1 f 1 du 



iy f 1 ----- 



m=0 Je 



-■m\r v 



sin ( ;-.,.( '"y) (r 2 e mXT ?y) 



1 f du . 

— / — sin(riu) sm(r 2 u) 
A Jo u 



-1m 

2A 



r\ + r 2 



n - r 2 



(A15) 



The fact that this integral diverges for n = r 2 is an artifact of choosing the same cut-off for both encounters. When 
the cutoffs differ the resulting expression is finite and of order 1/A. (Certainly, there are no factors of t p .) For the 
summation with I > m we find 



oo oo „i j 

A ^ f-' _L-xt p y 



m=Q k=l Je 



■ sin (Tie 



-. ' n \ J 

Kr *y) sin (r 2 e- im+k)XT "y) =-^ / — sin(rm) sin (r 2 e- kXr 

k=i ^° u 







(A16) 



Note that while the argument in the first sine is typically large this may not be the case for the second sine. Let 
us denote by fco the point where the exponential balances r 2 , that is, r 2 e~ k ° XTp — 1. Typically fco is not an integer 
and r 2 e~[ fc °l Arp 1 while r 2 e -( -[ fc °l +1 - )Arp <C 1. For k < fco both sines are rapidly oscillating and the integral can be 
approximated by integrating up to infinity. This result in the series 



1 



[ fc o] .oo 



fc=l 



du 



[fco] 



aE / — sin(riM)sin(r 2 e kXT ?u) ~ — In 



2A 



k=l 



r\ + r 2 e 



-k\r v 



T\ - r 2 e 



-k\T 



A 6 



-\t„ 



(A17) 



For k > fco we can expand the second sine since its argument is much smaller than unity everywhere in the integration 
range. As a result one finds 

e -([fco] + l))Ar p 



- y 

A f-' 



r 2 e 



„-/cAt 



fc=[fc ]+i 



sv[\(riu)du ~ — (1 — cosri) 
A 



1 



(A18) 



This part of the sum is much smaller than the other con- 



tributions. The case where m > I is treated in the same 
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FIG. 8: (Color online) In the presence of time-reversal sym- 
metry one also has to consider interference of time-reversed 
trajectories. The figure shows the configuration of type A 
after time-reversal of the trajectories 73 and 74. 

way. Thus, we can conclude that the series of integrals 
in Eq. (|AL4[1 is of typical size r p /A. 

Combining the integrals J\ and J 2 we find that J 2 — 
•h,\Ji,2 is of order r p /A. Since it is in this combination 
that these integrals enter into the Eq. (|A4[) for J, we 
conclude that J is of order td / A, as advertised. 

APPENDIX B: TIME REVERSAL SYMMETRY 

The calculation of the pumping current variance per- 
formed in Sec. IHII did not include contributions from tra- 
jectory configurations which depend on the presence or 
absence of time reversal symmetry. Indeed, all stretches 
of interfering trajectories in Fig. [3] are traversed in the 
same direction, so that any accumulated phase difference 
does not depend on an applied magnetic field. For con- 
ductance fluctuations, there is an equal contribution to 
the variance of the conductance in which one pair of tra- 
jectories is replaced by the time-reversed trajectories, so 
that the interfering stretches are traversed in opposite 
directions, as, e.g., in Fig. [8] (In the language of dia- 
grammatic perturbation theory, such contributions are 
referred to as 'Cooperons'.) Such contributions depend 



on the presence or absence of time-reversal symmetry and 
are responsible for the factor- two suppression of var G 
upon application of an external magnetic field.— ^ For a 
quantum pump, however, there is no contribution from 
trajectories that traverse the interference region in oppo- 
site directions, so that (I 2 ) is independent of the presence 
or absence of time- reversal symmetry. This is well known 
in terms of random matrix theory 11 ! 23 A semiclassical 
derivation of this fact is given below. 

The contribution of type A to the correlator (n i II ii ) 
is shown in Fig. [5J but with the trajectories 73 and 
74 replaced by their time-reversed. The main difference 
between this configuration of trajectories and configura- 
tions which do not depend on time reversal symmetry is 
that the trajectories in the central loop propagate in op- 
posite directions. We'll now show that the contribution 
of such a configuration of classical trajectories to (ll L Il n } 
is zero, even in the absence of an applied magnetic field. 

We first consider the contribution of the type "Al", 
in which both encounter regions fully reside inside the 
quantum dot. For this contribution, the reversal of the 
direction of two of the trajectories has no effect, and one 
finds that this contribution is given by Eq. (I20[) of Sec. 
HID For the configuration depicted in Fig. [8] both en- 
counters can touch the leads, resulting in three additional 
contributions to (H L H R ). Each of the configurations in 
which only one encounter touches the lead openings has 
a contribution given by Eq. (f22|) of Sec. IIIII (up to the 
interchange i cn c,i t cnc ^ 2 in one of these). The last con- 
tribution, from trajectory configurations in which both 
encounters touch the leads, was not considered in Sec. 
IIIII In such a configuration, 71 and 72 enter through 
the right contact and 73 and 74 to enter through the left 
contact. The calculation of this contribution is similar 
to that of the contribution "A2" in Sec. IIIII but with in- 
tegrations over two encounter times instead of only one. 
Labeling this configuration "A4" , we find 



(U L U R ) Ai = 32Nl % 27 » CiC 2 r 2 \ 2 I dx 1 dx 2 cos(rx 1 )cos(rx 2 ) (Bl) 



>A4 jV2 

x { 5 r D (l - x\' X ^) (l - x\ /X ^) - 2t cnc , 2 , 2 /ATD (l - 4 /X/T °) ~ 2 W ^ A - (l - x\' X ^) } , 

where t eac j = (1/A) \i\{\/xj), j — 1,2. Combining all contributions we then find that the additional contribution to 
the correlator (Yl L H R ) in the presence of time-reversal symmetry is 

(n L n«)7 S = ^j^^ CiC 2 r 2 \ 2 J' d Xl dx 2 cos(r Xl ) cos(ra 2 ) (2x\ /Xt » + 2x\ /Xt ° - h) . (B2) 



Upon integration over X\ and £2, all three terms between contributions of type B proceeds along in exactly the 
brackets yield fast oscillating functions of r — c 2 /h, which same way. We have verified that the resulting contribu- 
are neglected in the classical limit. The calculation of tion can also be neglected. 
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